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Abstract 

This  report  covers  technical  progress  on  the  Wide-field  Imaging  System  and  Rapid 
Direction  of  Optical  Zoom  (WOZ)  contract  for  the  period  16  Dec  2010  -  15  Mar  2011.  A  report 
covering  financial  status  for  this  period  is  provided  separately.  During  this  period,  the 
multiphysics  design  tool  has  been  further  developed  and  is  nearing  completion.  A  substantial 
amount  of  testing  on  the  initial  materials  samples  has  revealed  several  characteristics  which  are 
not  yet  well  understood.  The  set  of  material  characterization  samples  has  been  fabricated  by 
NeXolve  and  will  undergo  additional  processing  in  the  coming  quarter. 
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Summary 

In  June  2010,  Digital  Fusion  Solutions  (DFS,  a  wholly-owned  subsidiary  of  Kratos 
Defense  and  Security  Solutions,  Inc.)  was  awarded  a  contract  for  the  development  of  optical 
quality  flexible  membrane  mirrors.  This  report  covers  technical  progress  on  the  Wide-field 
Imaging  System  and  Rapid  Direction  of  Optical  Zoom  (WOZ)  contract  for  the  period  16  Dec 
2010  -  15  Mar  2011.  A  report  covering  financial  status  for  this  period  is  provided  separately. 
During  this  period,  the  multiphysics  design  tool  has  been  further  developed  and  is  nearing 
completion.  A  substantial  amount  of  testing  on  the  initial  materials  samples  has  revealed  several 
characteristics  which  are  not  yet  well  understood.  The  set  of  material  characterization  samples 
has  been  fabricated  by  NeXolve  and  will  undergo  additional  processing  in  the  coming  quarter. 


Introduction 

DFS  has  teamed  with  Advanced  Optical  Systems,  Inc  (AOS)  and  NeXolve,  Inc.  to 
develop  a  breadboard  system  for  demonstrating  non-mechanical  zoom  using  flexible  thin  films. 
The  project  consists  of  three  major  task  area,  which  are  Material  Characterization,  Diagnostic 
System  Development,  and  Imaging  System  Development.  DFS  is  responsible  for  leading  the 
modeling  effort,  AOS  is  responsible  for  leading  the  optical  design  effort,  and  NeXolve  is 
responsible  for  fabrication  of  the  films.  In  addition,  each  member  of  the  team  supports  the  others 
in  their  areas  of  responsibility. 

Material  Characterization  phase  will  involve  two  components:  development  of  modeling 
tools  and  measurement  of  material  properties  for  use  in  the  tools.  This  phase  is  currently 
ongoing  and  will  be  discussed  thoroughly  in  the  next  section.  The  Diagnostic  System 
Development  phase  will  apply  the  tools  and  knowledge  developed  in  the  first  phase  to  build 
algorithms  for  modeling  and  testing  one-dimensional  films  with  multiple  actuators.  These 
models  will  then  be  compared  to  experimental  results  and  modified  as  needed.  The  Imaging 
System  Development  phase  will  conclude  the  effort  by  building  a  sensor  system  with  variable 
zoom  based  on  optical  quality  PVDF  films. 

Methods,  Assumptions,  and  Procedures 

The  modeling  tools  are  based  on  interaction  between  three  commercial  software 
packages:  SolidWorks,  COMSOL  Multiphysics,  and  ZEMAX  optical  design.  The  multiphysics 
design  tool  is  nearing  completion.  We  have  demonstrated  the  ability  to  create  a  model  in 
SolidWorks,  import  the  model  into  COMSOL,  and  then  interactively  update  the  model  with 
either  package.  COMSOL  can  then  perform  the  electrostatic  and  mechanical  modeling  to 
calculate  the  deformation  resulting  from  the  applied  voltages.  Linally,  the  deformed  surface  can 
be  exported  to  ZEMAX  via  MatLab.  Lrom  ZEMAX,  various  analyses  can  be  conducted  to 
determine  important  parameters  such  as  focal  point,  aberrations,  and  wavefront  distortion. 
Defining  the  specific  analytical  results  to  extract  from  ZEMAX  to  support  the  optimization 
remains  to  be  determined. 

Analytical  and  Numerical  Model 

An  analytical  model  was  developed  to  determine  the  relationship  between  the  curvature 
of  a  clamped  cantilever  beam  and  the  voltage  applied.  A  model  for  a  bimorph  piezoelectric 
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Structure  from  the  literature  was  used  to  compare  the  results  with  the  numerical  results  from 


COMSOL.  From  Smits/  an  expression  for  the  displacement,  5,  is  given  by 
^  L^F  3^3  iLV 

S  = - ^  + 

2Ewh^  Sh^ 


(1) 


where  L  is  the  length  of  the  film,  w  is  the  width  of  the  film,  h  is  the  thickness  of  each  layer  of  the 
bimorph,  E  is  Young’s  modulus,  ^31  is  the  piezoelectric  coefficient  relating  the  strain  in  the  x 
direction  to  the  field  in  the  z  (thickness)  direction,  F  is  the  applied  force,  and  V  is  the  applied 
voltage. 

Table  1.  Comparison  of  numerical  and  analytical  results  for  bimorph  cantilever  with  applied  force  and 


voltage 


Inputs 

Outputs 

Displacement  (mm) 

Voltage  (V) 

Force  (N) 

Numerical 

Analytical 

20000 

0 

32.812 

32.812 

0 

-0.00152 

-32.716 

-32.803 

2000 

-0.00152 

0.097 

0.010 

Again,  following  Smits,  the  charge,  Q,  generated  on  each  surface  from  an  applied  force  is  given 
by 


<2- 


3d^,L^F 


(2) 


Substituting  S  = 


L^F  3^3  iLV 

2Ewh^  Sh^ 


with  the  voltage  set  to  zero  yields  charge  as  a  function  of  the  deflection: 


3d^^EwhS 

4L 


(Eq.  1 


(3) 


The  comparison  of  this  expression  with  the  COMSOL  results  is  given  in  Table  2  .  There  appears 
to  be  a  factor  of  2  difference  in  the  results.  This  difference  is  an  area  for  future  investigation, 
and  is  suspected  to  be  an  error  in  Reference  1 .  The  value  of  displacement  used  was  equal  to  the 
displacement  resulting  from  the  application  of  2  kV  with  zero  applied  force. 


Table  2.  Comparision  of  numerical  and  analytical  results  for  displacement 


Inputs 

Outputs 

Displacement 

(mm) 

Numerical 
Charge  (nC) 

Analytical 
Charge  (nC) 

32.812 

50.2 

25.0 

Given  the  generally  good  agreement  between  the  bimorph  analytical  expressions  and  the 
numerical  model,  we  compared  the  results  for  an  analytical  expression  of  displacement  resulting 
from  an  applied  voltage  to  the  corresponding  numerical  results.  To  our  knowledge,  this  is  the 
first  derivation  of  this  relationship.  Modeling  the  piezoelectric  constant  as  ^31  +  zd^\  in  a  single 
layer  unimorph,  we  found 

do, 

5  =  ^ — E.  (4) 

2h 
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Several  variations  of  voltage  and  film  thickness  were  tested,  and  all  results  agreed  to  at  least  five 
significant  figures.  As  an  example,  using  a  70  pm  thick  film  with  d-i\  =  5x10'  C/J  (a  value 
typical  of  what  we  expect  based  on  d^i  =  28x10'  C/N)  yielded  a  deformation  of  54.7  mm  at  the 
end  of  an  87.5  mm  film.  However,  this  simplified  model  does  not  account  for  the  internal  field 
that  results  from  the  strain  caused  by  the  bending  of  the  film.  When  the  numerical  model 
includes  this  effect,  the  deformation  is  reduced  by  approximately  14%.  This  effect  clearly 
warrants  further  investigation  and  we  will  attempt  to  include  it  in  the  analytical  model. 

The  multiphysics  design  tool  is  nearing  completion.  Previously,  several  models  of  one 
dimensional  film  with  varying  electrode  spacing  were  simulated.  By  examining  the  curvature, 
the  derivative  of  the  deformation,  along  the  surface  of  the  electrode,  the  optimal  spacing 
between  electrodes  can  be  determined.  The  result  was  that  the  electrodes  need  to  be  spaced  1mm 
or  less  to  maintain  a  continuous  curvature  without  any  drastic  deviations.  For  the  2-D  model,  the 
electrodes  will  be  in  a  pattern  consisting  of  rings  and  arms  to  set  the  desired  curvature.  The 
optimization  of  the  electrodes  requires  an  algorithm  developed  to  determine  the  locations  of  the 
electrodes  based  on  inputs  such  as  the  number  of  rings  and  arms,  electrode  spacing,  and  size  of 
the  outside.  MATLAB  was  used  to  develop  the  algorithm  since  it  is  the  central  hub  of  the 
integrated  software  model.  Figure  1  is  the  result  from  a  single  run  of  the  algorithm. 


270 


Figure  1.  Example  electrode  pattern  from  script  within  the  multiphysics  tool 

Experimental  Testing 

Testing  has  taken  two  separate  directions.  The  first  area  is  exploration  of  the  deformation 
response  of  the  film  to  an  applied  voltage,  and  the  approach  was  described  in  the  previous 
quarterly  report.  The  films  have  been  subjected  to  a  broad  series  of  tests  described  below.  The 
second  area  is  exploration  of  the  charge  generation  response  of  the  film  to  an  applied  load.  Both 
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tensile  and  bending  loads  have  been  applied  and  response  measured.  However,  the  bending 
loads  have  not  been  as  well  controlled  so  the  results  are  not  as  repeatable. 

Applied  Voltage 

During  this  period  the  experimental  setup  included  a  digital  to  analog  convertor  (DAC), 
which  enabled  computer  control  of  the  applied  voltage  as  well  as  measurement  of  the 
deformation  as  previously  described.  One  of  the  first  tests  was  application  of  a  0  -  2  kV  square 
wave,  which  is  shown  in  Figure  2  (top).  The  pronounced  ringing  was  due  to  the  mechanical 
vibration  of  the  excess  portion  of  the  film  below  the  electrode.  After  trimming  the  film  to 
remove  the  excess,  the  ringing  was  virtually  eliminated  Figure  2  (bottom). 


Figure  2.  Tip  response  to  alternating  0  and  2kV  voltage,  10.6  cm  length  (above)  and  3.8  cm  length  (below) 

From  this  result,  it  can  be  seen  that  there  is  a  significant  amount  of  drift  in  the  film.  The 
remaining  data  taken  during  this  test  (which  lasted  almost  900  sec)  show  the  drift  continues  until 
the  deformation  at  0  V  is  almost  0.5  mm.  A  similar  pattern  was  observed  when  a  constant  3  kV 
was  applied  to  the  film.  Repeated  tests  have  produced  somewhat  inconclusive  and  inconsistent 
results.  These  results  include  what  appears  to  be  a  spontaneous  reversal  of  the  poling  direction 
of  the  film,  significant  change  in  deformation  in  the  absence  of  applied  voltage,  and  changes  in 
the  responsivity  of  the  film  to  applied  voltage.  We  have  engaged  in  a  systematic  study  to 
understand  the  cause  of  these  results  and  correlate  the  response  with  the  experimental  conditions. 
Initial  indications  include  temperature  of  the  film  may  be  a  factor.  Measured  temperatures 
reached  46  °C.  However,  published  results  indicate  thermal  effects  on  the  piezoelectric 
coefficients  are  fairly  minor  until  temperatures  of  80-140  °C  (depending  on  the  polymer)  are 
reached.^  Testing  using  LED  illumination  rather  than  incandescent  illumination  is  being 
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performed,  but  thermal  effects  will  need  to  be  revisited.  It  also  appears  that  polymer  creep  may 
play  a  role  in  the  behavior. 


Applied  Load 


The  applied  load  testing  has  used  an  ORTEC  439  Digital  Current  Integrator  to  measure 
the  charge  generated  by  applying  various  loads  to  the  films.  Two  configurations  have  been  used 
in  which  the  applied  loads  have  been  in  the  plane  of  the  film  (tensile  test)  and  perpendicular  to 
the  surface  of  the  film  (bending  test).  The  tensile  testing  has  been  performed  using  the  NeXolve 
Instron  tensile  tester.  An  example  is  shown  in  Figure  3.  As  the  jaws  holding  the  ends  of  the  film 
separate  vertically,  charge  is  generated  on  the  two  electrodes  by  the  piezoelectric  response  of  the 
film.  The  charge  is  collected  by  the  ORTEC  439,  which  outputs  a  digital  pulse  for  each  100  pC 
of  charge.  These  pulses  are  counted  and  recorded  as  a  function  of  time.  Simultaneously,  the 
Instron  is  recording  the  applied  load  and  extension  of  the  film  as  a  function  of  time.  By 
synchronizing  the  times  (the  end  of  the  test  was  the  sharpest  event),  the  charge  vs.  load  curve 
was  determined. 


Figure  3.  Film  loaded  in  Instron  tensile  tester  with  connections  to  the  current  integrator 
Following  the  constituent  equations  given  in  Smits, 

1  11  1  31  3  ^ 

E3  and  D3  are  the  electric  and  displacement  fields,  respectively,  in  the  z-direction  (thickness  of 
the  film),  and  5  n  is  the  compliance  in  the  x-direction  (long  direction,  which  is  the  direction  of 
the  load).  Stress  (7)  and  strain  (S)  are  given  by 
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Ax  F 

s,=  —  r,  =  — 

L/  hw^ 


(6) 


with  symbols  Lf,  length  of  film;  w/,  width  of  film;  h  thickness  of  film;  Lg,  length  of  electrode;  We, 
width  of  electrode;  F,  applied  load;  and  Ax,  extension  of  film  under  load. 

With  the  current  integrator  connected,  the  electrodes  are  effectively  connected  due  its 
very  small  input  impedance,  so  £3  =  0,  and 


L  V 

Ax  =  ^F 

hWf 


hw, 


or 


(7) 


where  nix-F  is  the  slope  of  the  extension-load  curve.  The  compliance  is  equal  to  the  reciprocal  of 
Young’s  modulus  (£),  which  is  often  referenced.  An  example  set  of  data  is  shown  in  Figure  4. 
As  expected,  the  result  is  very  linear.  The  elastic  modulus  determined  from  this  graph  is  1.52 
GPa,  which  is  reasonably  consistent  with  average  literature  values  of  PVDF  of  1.65  GPa. 


The  bound  surface  charge  density  developed  as  a  result  of  the  strain  is  oz,  =  P,  where  P  is 
the  polarization."^  By  definition,  D  =  P  +£oE,  so  Oh  =  dsiF  =  Q/Ag  =  QILgWg.  Here  Q  =  kC,  where 
k  is  the  current  integrator  setting  (100  pC/count)  and  C  is  the  number  of  counts.  Thus, 

Q  =  d^\T\  Ag=  d^i{LgWg)l{LfWf)  (8) 


and 


C  -  d 


31 


khWf 


or 


— 


hw  ^ 
kLw„ 


-m, 


C-F  ’ 


(9) 


where  mc-F  is  the  slope  of  the  counts  vs.  load  curve.  The  pulse  count  vs.  load  graph  is  also  very 
linear,  as  shown  in  Figure  4.  The  calculated  value  of  d^i  is  12.5  pC/N  (or  pm/V).  This  is 
somewhat  lower  than  the  values  seen  in  the  literature  of  15-28  pC/N,  but  very  reasonable. 
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Figure  5.  Current  integrator  counts  vs.  applied  tensile  load  for  the  same  test  as  Figure  4 

The  first  round  of  testing  was  primarily  to  explore  the  test  procedure.  Based  on  this 
exploratory  testing,  we  have  identified  the  appropriate  test  conditions,  such  as  strain  rate.  This 
testing  occurred  near  the  end  of  the  quarter,  and  additional  tests  are  scheduled  for  the  first  week 
of  next  quarter. 

The  second  set  of  tests  with  applied  loads  involved  applying  a  lateral  load,  which  resulted 
in  the  film  bending.  This  was  preliminary  testing  to  explore  test  procedures  and  shown  in  Figure 
7.  As  before,  testing  was  started  using  the  bimorph  sample  fabricated  from  two  films  bonded 
together  with  poling  directions  anti -parallel.  From  Eq.  1,  a  linear  relationship  is  expected 
between  the  displacement  and  the  charge  generated.  However,  Figure  8  shows  that  there  is  a 
quadratic  component  to  the  charge  production.  This  could  be  due  to  the  presence  of  the  gradient 
in  the  piezoelectric  coefficient,  since  we  have  not  yet  developed  an  analytical  model  for  the 
charge  generation  including  the  gradient,  as  well  as  other  effects  not  accounted  for.  However, 
ignoring  the  quadratic  component  for  present,  a  linear  fit  using  the  measured  value  of  E  =  1.6 
GPa  yields  a  d^\  value  of  11  pm/V,  which  is  in  agreement  with  the  value  found  from  the  tensile 
testing.  Using  the  slopes  from  the  upper  half  of  the  curve  and  the  lower  half  curve,  dsivalues  in 
the  range  8-15  pm/V  are  obtained. 

Dielectric  Constant 

Although  the  dielectric  constant  is  not  used  in  the  analysis  of  the  experimental  data,  it  is 
used  in  the  COMSOL  modeling.  Therefore,  a  quick  measurement  of  the  dielectric  constant  was 
performed.  A  circuit  was  built  as  shown  below.  The  AC  voltage  source  was  the  NI 9264  analog 
output.  The  source  was  configured  to  output  0  -  10  V  square  waves.  The  voltage  measurement 
was  done  using  the  Agilent  oscilloscope.  The  resistor  was  the  49.9  kohm  precision  resistor.  The 
capacitor  was  the  69  pm  film  sample  E. 
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Figure  6.Schematic  diagram  for  capacitance  measurement 

The  oscilloscope  recorded  the  decay  of  the  voltage  as  the  voltage  source  transitioned 
from  10  V  to  0  V.  After  allowing  transient  from  analog  output  to  decay  for  20  ps,  voltage  vs 
time  is  shown  in  the  figure  below.  The  exponential  decay  provides  excellent  fit  to  the  data. 


Capacitance  measurement 


time  (s) 


From  r  =  R/C,  the  capacitance  can  be  determined  to  be  731  nF. 
capacitance  of  a  parallel  plate  capacitor, 


c  =  ^,^0 :7 
a 


Then  using  the 


(10) 


and  the  dimensions  of  the  electrodes  (L=  52  mm,  w=  11  mm,  d=69  pm),  the  dielectric  constant 
£';^10.0.  This  is  in  reasonable  agreement  with  the  value  found  in  the  literature. 
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Figure  7.  Test  setup  measuring  response  to  lateral  loads 


Figure  8.  Charge  generated  by  bending  bimorph  strip 


Material  Crystallinity 

Films  manufactured  during  the  previous  two  quarters  were  examined  further  by  NeXolve 
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using  differential  scanning  calorimetry  (DSC)  to  quantify  the  relative  amounts  of  beta  and  alpha 
crystal  phases  present  in  the  film  samples.  The  beta  crystal  phase  is  piezoelectrically  active 
while  the  alpha  phase  is  not.  DSC  can  determine  the  relative  amounts  of  each  crystal  phase, 
although  NeXolve  is  not  currently  set  up  to  determine  the  absolute  material  fraction  of  each 
phase.  The  proportions  of  each  crystal  phase  were  measured  and  can  be  seen  in  Table  3. 


Table  3:  Relative  amounts  of  Beta  and  Alpha  Crystal  phases  in  PVDF  copolymer  films. 


Sample 

Film 

Thickness 

Beta 

Crystal 

Alpha 

Crystal 

Total 

Crystallinity 

Beta 

Crystal 

Alpha 

Crystal 

Identification 

um 

J/g 

J/g 

J/g 

% 

% 

081810-1010 

54 

31.87 

62.86 

94.73 

33.64% 

66.36% 

081810-1012 

55 

64.41 

27.95 

92.36 

69.74% 

30.26% 

081810-1014 

55 

33.99 

61.71 

95.70 

35.52% 

64.48% 

BF100908.004 

41 

49.34 

6.295 

55.63 

88.69% 

11.31% 

BF100908.003 

75 

48.95 

8.449 

57.39 

85.28% 

14.72% 

BF100908.005 

69 

49.05 

6.425 

55.47 

88.42% 

11.58% 

aftertesting- 

BF100908.005-A 

69 

50.05 

6.551 

56.60 

88.43% 

11.57% 

aftertesting- 

BF100908.005-E 

69 

50.12 

7.212 

57.33 

87.42% 

12.58% 

The  first  three  samples  (green)  are  examples  of  films  from  other  programs  that  have  the 
incorrect  amount  of  beta  and  alpha  crystal  phases  needed  to  be  piezoelectrically  active,  while  the 
other  samples  are  films  produced  for  this  program  that  exhibit  the  proper  crystal  phases.  The 
samples  in  orange  are  from  the  same  film.  Sample  BF100908.005  was  tested  before  poling,  and 
after  poling  and  electrical/mechanical  testing  was  performed.  These  films  did  not  exhibit  a 
significant  change  in  crystallinity  before  poling  and  after  testing. 


Results  and  Discussion 

The  multiphysics  design  tool  is  integrating  well  and  will  enable  an  empirical  approach  to 
understanding  the  behavior  of  the  films.  The  next  step  in  the  tool  development  is  determination 
of  the  parameters  to  extract  from  the  optical  analysis  in  order  to  compute  the  voltages  needed 
achieve  the  desired  shape.  This  task  will  be  closely  tied  with  development  of  the  diagnostic 
system  optical  design,  which  will  begin  early  in  the  next  quarter.  The  need  for  an  analytical 
model  for  charge  generated  given  an  applied  load  or  deformation  is  needed  for  understanding 
behavior  of  the  film  and  under  development. 

On  the  experimental  side,  the  biggest  challenge  is  developing  an  understanding  of  the 
long  term  behavior  of  the  films  under  applied  voltages.  In  addition,  the  effect  of  relatively  low 
heat  levels  is  not  consistent  with  our  expectations  based  on  results  from  the  literature.  We  have 
begun  a  systematic  study  to  determine  the  conditions  and  causes  of  the  anomalous  behavior.  The 
results  obtained  by  tensile  testing  have  been  consistent  with  values  from  the  literature  as  well  as 
values  from  other  testing  performed  on  this  program.  The  next  step  will  be  assessing  a  variety  of 
sample  under  better  defined  conditions.  The  bending  tests  will  be  performed  using  the  unimorph 
films  with  a  motorized  translation  stage  to  assess  the  impact  of  translation  speed  on  the 
measurements  and  to  improve  repeatability. 
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Conclusions 

The  multiphysics  design  tool  is  progressing  according  to  plan.  At  the  beginning  of  the 
quarter,  the  experimental  work  was  behind  schedule  but  is  catching  up.  However,  there  have 
been  unexpected  challenges  in  understanding  the  properties  of  these  materials.  The  focus  for  the 
next  quarter  will  be  on  completing  the  material  characterization  and  beginning  the  design  of  the 
diagnostic  optical  system. 
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